The Kondo effect in quantum dots at high voltage: Universality and scaling 
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We examine the properties of a dc-biased quantum dot in the Coulomb blockade regime. For voltages 
V large compared to the Kondo temperature Tk , the physics is governed by the scales V and 7, where 
7 ~ V / \v?(V /Tk) is the non-equilibrium decoherence rate induced by the volt age- driven current. 
Based on scaling arguments, self-consistent perturbation theory and perturbative renormalization 
group, we argue that due to the large 7, the system can be described by renormalized perturbation 
theory for \xi{V/Tk) 3> 1. However, in certain variants of the Kondo problem, two-channel Kondo 
physics is induced by a large voltage V . 
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In recent years, it became possible to observe the 
Kondo effect in quantum dots in the Coloumb blockade 
regime &-0|. These systems allow to investigate, how 
non-equilibrium induced by external currents and bias 
voltages influences the Kondo physics. Similarly, the ex- 
perimentally observed anomalies of the energy relaxation 
in strongly voltage-biased mesoscopic wires || have re- 
cently been shown |(| to be caused by scattering from 
magnetic impurities or two-level systems. 

In equilibrium, almost all properties of the Kondo 
effect are well understood, and the Kondo model to- 
gether with the methods used to solve it (e.g. renormal- 
ization group (RG), Bethe ansatz, conformal field the- 
ory, bosonization, density matrix RG, flow equations or 
slave particle techniques) have become one of the central 
paradigms in condensed matter theory. However, in non- 
equilibrium many of the above-mentioned methods fail, 
and despite the experimental and theoretical relevance 
and a substantial body of theoretical work ]7|-|l5[|, sev- 
eral even qualitative questions about the Kondo effect in 
non-equilibrium have remained controversial. Recently, 
Coleman et al. |l4|] claimed that the Kondo model at high 
voltages V ^S> Tk cannot be described by (renormalized) 
perturbation theory (PT) but is characterized by a new 
two-channel Kondo fixed point (see also |l3) ). By con- 
trast, Kaminski et al. || argue that the non-equilibrium 
decoherence rate 7 destroys the Kondo effect. Wc will 
show in the following that the Kondo effect is indeed de- 
stroyed in the case of the usual Anderson model, but for 
certain variants of the Kondo model, where the current 
at high bias is suppressed, the scenario proposed in |l4| ] 
appears to be recovered. 

We model the quantum dot using the Anderson model 



H A = H +e d ^2dld a + ^2{t a cl kry d a + h.c.) 
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where H = J^aka ^ak^ aka c ak<y is the Hamiltonian of the 
electrons in the left and right leads, a — L, R, charac- 
terized by a dc bias voltage V, e L / Rk — e k ± V/2, re- 



spectively. We will consider only symmetrical dots with 
tunneling matrix elements tj, = t R = t. The negative Ed 
with \ed\ ^> r = 2irNot 2 , where iVo is the electron density 
of states in the leads, and the large Coloumb repulsion 
U — ► 00 enforce the number of electrons rid — J2a ^t^o- 
in the dot level to be approximately 1. 

In this regime, the local degree of freedom of the quan- 
tum dot is a spin S = \T1 I<J d\ r a aa ,d <y i , where a is 
the vector of Pauli matrices, and the low-energy proper- 
ties of Ha are well described by the two-lead Kondo (or 
Coqblin-Schrieffer) model, 

H K =H + V Q J>[ 0ff + 4j(C£0<r + CflOa) + (2) 

(T 

where cl/rq<j — J2k c L/Rka- For an Anderson model 
with symmetrical coupling to the leads, one obtains 
Jl = Jr. = Jlr. = Jrl = 4V - 2t 2 /e d = J. For 
sufficiently small J, the potential scattering term Vo can 
be neglected and, as will be seen, the equilibrium and 
the non-equilibrium physics of the Kondo model is com- 
pletely universal, characterized by a single scale, the 
Kondo temperature T K = D^J N Je 1/{2NaJ) , where D 
is a high-energy cutoff. The precise formula for the pre- 
factor of Tk depends on details of the model. However, 
for Tk ,T,V <C D relevant physical quantities like the 
conductance G are universal, G = G(V/Tk,T/Tk) and 
do not depend on details of the original Hamiltonian. 

In the first part of the paper, we investigate in detail 
the Anderson model in the Kondo regime at high voltages 
using the so-called non-crossing approximation (NCA). 
In the second part we will use the insight gained from this 
analysis to study a heuristic version of poor man's scaling 
in non-equilibrium for a Kondo model with Jlr < Jl/r- 

To derive NCA, one first rewrites Ha in the limit U — > 
00 using a so-called pseudo-fermion f a and a spin-less 
slave boson b with d a — W / CT , subject to the constraint 
Q = J2a flfv + b'b = 1. The Anderson model then takes 
the form H A - H + e d tfb + £ Q jV a cl aQ tff a + h.c). 
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In this language, the NCA is just the lowest-order self- 
consistent PT in t a , where the constraint Q = 1 is taken 
into account exactly. While the NCA fails to describe the 
low-energy Fermi liquid fixed point in the Kondo regime 
correctly |l(|, it gives reliable results (with errors of the 
order of 10%) in equilibrium for temperatures down to a 
fraction of Tk- As a self-consistent and conserving ap- 
proximation, it also displays the correct scaling behavior 
and reproduces the relevant energy scales. 

While the NCA equations in non-equilibrium have 
been solved by many groups |l5fl , we are not aware of any 
careful analysis of the relevant scales at high bias voltage, 
which is central for a qualitative understanding of the 
non-equilibrium Kondo effect. Generally, the NCA equa- 
tions have to be solved numerically; however, in the limit 
of extremely high voltage, In V/Tk ^S> 1 (but ^<B),an 
analytical solution is possible: the problem is in the weak 
coupling regime. Finite V induces an inelastic spin relax- 
ation or decoherence rate. Since in NCA the spin density 
is just a convolution of the pseudo-fermion propagators, 
this rate is given by 2Im£/(0) = 27, with £/ the pseudo- 
fermion self-energy. We start by calculating the retarded 
self-energy ££(w) of the boson, using the fact that (as 
shown below) the spectral function of the pseudo-fermion 
is a sharp peak of width 7 <C V . Throughout we consider 
the low temperature limit, T = 0, and obtain Im££(u;) 
-7rJiV | £d |(/2(-cj) + /^(-ct))), where f R/L are the Fermi 
functions in the left and right leads, broadened by 7. The 
step- function in ImEj(w) leads to logarithmic contribu- 
tions to ReEJ(oj), cut off by 7 and the band width D. Us- 
ing relations like 1 - 2N J\n[D/\w\] = 2N J\n[\u\/T K ] 
one obtains for the real part of the boson propagator 
G[(w), for lnV/T K > 1, 



No J? b ca (lu) = 2N t 2 BeG b (w) 
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This combination plays the role of an effective (frequency 
dependent) exchange coupling J c g. Remarkably, the per- 
turbative expression Eq. (0) would develop a pole close 
to lu = ±V/2 if 7 < T* = y/T'x + {V/2) 2 -V/2 « T 2 /V. 
The breakdown scale T* of PT has also been discussed in 
[4th Ref. of §] and @. It indicates that Eq. (§) is only 
valid for T* < 7 < V. Indeed, this criterion is fulfilled 
(see Fig. [I]), as one finds within NCA 
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For the conductance in units of the conductance quantum 
G = 2e 2 /(2irh) we obtain for lnV/T K > 1 
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FIG. 1. Non-equilibrium decoherence rate 7 = ImE/(0) 
calculated within NCA compared to the strong coupling scale 
T* = ^Tf+WJW ~ V/ 2 (dashed line). For 7 > T* one 
stays in the weak coupling regime. The symbols correspond to 
3 different values of Tk, T k /D = 9- 1(T 6 , 8- 1(T 5 , 5- 10" 4 . In- 
set: conductance G in units of Go = 2e 2 /(2nh). Long-dashed 
and solid lines: asymptotic analytical results, Eqs. (^), (^|), in 
leading and next-to-leading order, respectively. 

Numerical results for smaller voltages down to V < Tk 
are shown in Fig. [I] and display universal behavior over 
the complete range of voltages and over several orders of 
magnitude in Tk. Despite the fact that for high voltages, 
In V/Tk ^> 1, one stays in the weak coupling regime, the 
prefactors of ■f NCA and G NCA are not exact, since the 
NCA for the Anderson model treats the potential scat- 
tering Vo and the Kondo coupling J incorrectly on equal 
footing. It is not difficult to obtain the correct asymp- 
totic prefactors || by calculating 7 and G in leading or- 
der PT in J for the Kondo model Eq. §) (with V = 0) 
and by replacing J by 1/ (2 In V/Tk)- This corrects the 
leading term of the NCA results Eqs. (Q), @ by a pref- 
actor 3/4. It is, however, important to stress, that the 
asymptotic result in the limit In V/ Tk — > 00 is almost 
useless as, due to the logarithmic dependence (Eq. (|j)), 
sub-leading corrections are very large (e.g. still 10% for 
V/Tk = 10 6 ). 

In the limit of large V, the scale 7 influences quanti- 
ties like the conductance, where all electrons in an energy 
window V contribute, only slightly. The situation is dif- 
ferent for the local spin susceptibility on the quantum 
dot, x or the spectral function Ad(uj) of the elec- 

tron on the quantum dot. Ad{ui) calculated numerically 
within NCA is shown in Fig. |[ Like many groups before, 
we obtain two well defined peaks at voltages ±V/2. In 
the limit \uV/Tk — > 00 we find approximately A^ CA (u) 
w (7r 2 /r) [N J^ s ca (lo)] 2 1 with large but universal sub- 
leading corrections and a non-universal, (almost) con- 
stant potential scattering background of 0(r/e?). NCA 
incorrectly treats potential and spin flip scattering on 
equal footing and, thus, overestimates the asymmetry of 
the peaks w.r.t. u> —lu in the small J limit. This can 
be seen from an analysis of the Schrieffer- Wolff transfor- 
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FIG. 2. Spectral function A^ ca (ll>) for various voltages 
V, each calculated in NCA at two values of Tk (solid and 
dashed lines) differing by a factor of 10. The NCA strongly 
overestimates the anisotropy of the peaks. Inset: asymptotic 
behavior of A^ CA (u>) in the Kondo scaling limit. Squares: 
numerical NCA result; solid line: asymptotic expression for 
\n(V/T K ) -> oo. 

mation which shows that this asymmetry is non-universal 
and of O(NqVq). Since the antisymmetrical in u> contri- 
butions to Ad{oj) cancel in the integral for the conduc- 
tance, non-universal corrections to G are much smaller. 
Note that the logarithmic cusps of Ad(u>) have an ad- 
ditional, small rounding of 0(7) compared to Eq. (||), 
but for large voltage, the half width at half-maximum 
of the peaks is not given by 7 but by \p^fV /2 ps 
O.SV/\n(V/Tx) > 7 (see Fig. |). 

Our analysis of the Kondo model suggests, that qual- 
itatively different behavior can be expected, if the non- 
equilibrium relaxation rate 7 is sufficiently small, 7 < 
T K /V . Since a non-zero 7 requires finite current, e.g. 
within bare PT 7 cx JlrV, it is therefore interesting 
to study the Kondo model Eq. (||) for Jlr <C Jl,Jr, 
using ideas from perturbative RG. Such a model can- 
not be derived from a simple Anderson model but may 
arise in more complicated situations. Not much is known 
about how the concepts of fixed points and renormaliza- 
tion group can be applied to a non-equilibrium situation 
(see, however, pj]|). The problem is that in the pres- 
ence of a finite bias voltage, many physical quantities 
like the conductance are not determined by low-energy 
excitations even at T = 0, since all states with energies 
of order of the applied voltage V contribute. Therefore, 
a controlled perturbative RG must probably be formu- 
lated for the full frequency-dependent vertices in Keldish 
space. We will not try to develop such a method here but 
propose to use a heuristic version of poor-man's scaling 
adapted to the present situation. As usual, we inves- 
tigate, how coupling constants change, when the cutoff 
A of the theory is modified. As long as the cutoff is 
large compared to the voltage, we expect that the usual 
poor-man's scaling equations hold. For the model with 



N J L = N J R 
obtains [SO 



g d and N J LR = g LR and V < A one 
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These are the RG equations of a channel-asymmetric two- 
channel Kondo model, where the even and odd channels 
couple to the spin with coupling constants g e = gd + gLR 
and g = g d — gLR < <?e- Note that for the Anderson 
model, the odd channel decouples and g = 0. Two 
parameters, Tx and a, determine the physics of the 
channel-asymmetric two-channel Kondo model, 



Tj 
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(gd - gLR.){gd + gLR.) 
2gLR 



, (7) 



where Tx is defined by g e (Tx) = 1. The dimensionless 
number a is the natural parameter to characterize the 
channel anisotropy, since it is invariant under the per- 
turbative RG flow Eq. (||), i.e. a(A) = ao = const, for 
A > V. If higher orders of g are included in Eqs. (||), the 
prefactor of Tx, Eq. (0), changes and the definition of 
the RG invariant a has to be slightly adjusted (a dimen- 
sionless invariant characterizing the flow will exist even 
in higher orders). For the usual one-channel Kondo effect 
or the Anderson model Eq. (|l]), cvo = 0, while for ao — > 00 
the model is just the well-known channel-symmetric two- 
channel Kondo model. We will therefore investigate, how 
a will change for A < V in order to determine if the sys- 
tem flows towards a two-channel fixed point, that has 
been proposed by Wen and Coleman et al. [Q. For 
V Tx, Eq. (0) is valid down to A = V and we obtain 



9d(V) = \ 
9lr(V) = I 
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For A < V, the calculation of the RG flow is less obvious. 
Some of the logarithmically diverging vertex corrections 
of J are cut off by the voltage V, changing the RG flow 
to 



dgd 
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in complete agreement with the analysis of Coleman et 
al. |fl4|| . However, all remaining logarithmic contributions 
are cut off by the decoherence rate 7 as it is evident, e.g., 
from our analysis of NCA. Thus, Eqs. ( |Io[ ) are only valid 
for 7 < A < V. 

Since in the perturbative regime of the RG the bare 
coupling constant No Jlr is replaced by the renormalizcd 
one, gLR, we find 

J~Vg 2 LR (V) : (11) 
which is V/[2hx(V/T K )} 2 for a < l/\n(V/T K ) and 
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FIG. 3. 7/T* as a function of ao at V — WTk- Inset: a* as 
a function of ao for various voltages. For 7/T* <C 1, the sys- 
tem displays strong coupling behavior for T < T* . q* 3> 1 in- 
dicates that this regime is dominated by two-channel physics. 



Nevertheless, bare perturbation theory cannot be applied 
and even the leading order of renormalized perturbation 
theory does not give precise results in the experimentally 
accessible regime due to large sub-leading corrections. In 
variants of the Kondo model with Jlr < Jl, Jr> a large 
voltage can, however, induce a qualitatively new behavior 
reminiscent of two-channel Kondo physics. 
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man, C. Hooley, L. I. Glazman, W. Hofstetter, O. Par- 
collet and A. Zawadowski for valuable discussions, the 
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work (J.K.) and the Emmy Noether program (A.R.) and 
SFB 195 of the DFG for financial support. 



V/[4alln 4 {V/T K )} for a > l/\n[V/T K ]. (The precise 
prefactor is irrelevant for our discussion). If we assume 
for the moment that 7 is small, we find that g<i flows to 
strong coupling at a scale T* defined by gd(T*) = 1, 



T* w T K 



Tk 
V 



l/[l+Jaoh(7/Tjt]] 



(12) 



For ao — this scale coincides with the one intro- 
duced in |l4|] , where the effects of 7 have been ne- 
glected. The system will, however, only flow to strong 
coupling if 7 < T* , while it remains in the weak cou- 
pling regime for 7 ^> T*. For the usual Kondo- or An- 
derson model with ao = 0, 7 is always larger than T* 
for V > T K (as V/T K > ]nV/T K ), and we therefore 
conclude in contradiction to Ref. |Q that in the sym- 
metrical Kondo model there is no strong coupling regime 
for V 3> Tk- The situation is, however, different in the 
asymmetric model with ao > 1/2 (Fig. ||). Here the ratio 
7 /T* sa {V/T K )/{Aal\n 4 [V/T K ]) is small for V < V* « 
T K (4a 2 In 4 [4a 2 In 4 [4a 2 In 4 [...]]]), e.g. V* « 6 • lO 4 ^ 
for ao = 1. What is the nature of this strong coupling 
regime which is reached for Tk < V < V* and ao > 1/2? 
Insight into this question can be gained from a calcula- 
tion of a (A = T*), defined in Eq. (0). Note that in the 
regime 7 < A < V, a is not invariant under the RG flow 
Eqs. ©. We obtain 



a(T*) « In — I 1 + a In — ) w a In 2 — . (13) 
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Obviously, a is strongly enhanced by the voltage (e.g. 
ln 2 [10 3 ] ss 50). Since for a — > 00 the system maps 
to a two-channel Kondo problem, we conclude that for 
a > 1/2 and T K < V < V* the system will likely be 
dominated by two-channel physics over a large regime. 

In this paper, we have shown that the usual Kondo 
model at high voltages Tk < V < D is a weak coupling 
problem because relaxation processes allowed in non- 
equilibrium even at T = destroy the building up of the 
Kondo effect, as their decoherence rate 7 ^> T* w T K /V. 
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